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ABSTRACT 

There are several astrophysical situations where one needs to study the dynamics of 
magnetic flux in partially ionized turbulent plasmas. In a partially ionized plasma 
the magnetic induction is subjected to the ambipolar diffusion and the Hall effect in 
addition to the usual resistive dissipation. In this paper we initiate the study of the 
kinematic dynamo in a partially ionized turbulent plasma. The Hall effect arises from 
the treatment of the electrons and the ions as two separate fluids and the ambipolar 
diffusion due to the inclusion of neutrals as the third fluid. It is shown that these 
nonideal effects modify the so called a effect and the turbulent diffusion coefficient 
/3 in a rather substantial way. The Hall effect may enhance or quench the dynamo 
action altogether. The ambipolar diffusion brings in an a which depends on the mean 
magnetic field. The new correlations embodying the coupling of the charged fluids and 
the neutral fluid appear in a decisive manner. The turbulence is necessarily magneto- 
hydrodynamic with new spatial and time scales. The nature of the new correlations is 
demonstrated by taking the Alfvenic turbulence as an example. 
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1 INTRODUCTION 

The kinematic dynamo has revealed many an essential work- 
ings of an astrophysical dynamo for the generation of mag- 
netic field in objects varying from stars to molecular clouds 
to accretion disks. The kinematic dynamo (Parker 1955; 
! Steenbeck, Krause and Radler 1966; Moffatt 1970; Stix 
1972) is based on the possible generation of an electromo- 
tive force parallel to the mean magnetic field in a reflex- 
ionally asymmetric turbulence, the so called a effect. Here 
q is a measure of the net kinetic helicity. The correspond- 
ing turbulent diffusion coefficient j3 becomes a function of 
the mean turbulent kinetic energy. The scale separation is 
an integral part of the kinematic dynamo. A weakly ionized 
plasma is defined by the condition (Alfven and Falthammer 
1962) that the electron-neutral collision frequency v en ~ 
10- 15 n nV /8A B T/(7rm en ) is much larger than the electron- 
ion collision frequency u e i ~ 6 x lO'^mAZ 2 (KbT)~ 3/2 . 
This translates into the ionization fraction n c /n n < 5 x 
10 -n T 2 (Alfven and Falthammer 1962) where n's are the 
particle densities and T is the temperature in Kelvin. Al- 
though the ideal magnetohydrodynamics (MHD) is often 
used as a starting point of an astrophysical investigation, 
there are many systems with a rather low degree of ioniza- 
tion dominated by the charged particle-neutral collisions and 



the neutral particle dynamics. A major part of the solar pho- 
tosphere (Leake & Arber 2006; Krishan & Varghese 2007), 
the protoplanetary disks (Krishan & Yoshida 2006) and the 
molecular clouds (Brandenburg & Zweibel 1994) are some of 
the examples of weakly ionized astrophysical plasmas. The 
dynamo action in such a plasma would be affected by the 
multifluid interactions. The issue of possible disconnection 
between the sub-surface and the surface solar magnetic field, 
recently emphasized by Schiissler (2005), may have some 
bearing on the neglect of the neutral fluid-plasma coupling in 
the flux transport on the solar photosphere. Zweibel (1988) 
studied the dynamo process in a partially ionized plasma 
within a single fluid description. Including only the ambipo- 
lar diffusion, she determined the velocity, the density and 
the magnetic field fluctuations self consistently in the form 
of magnetohydrodynamic waves and thus went beyond the 
kinematic dynamo. We develop a three fluid framework for 
a kinematic dynamo including the Hall effect and the am- 
bipolar diffusion in section two. The a effect of the kinematic 
dynamo is formulated in section three. The new correlations 
arising due to the coupling amongst different fluids are un- 
derstood by taking the Alfvenic turbulence as an example 
and we end the paper with a section on discussion and con- 
clusion. 
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2 THREE-COMPONENT MAGNETOFLUID 

We begin with the three component partially ionized plasma 
consisting of electrons (e), ions (i) of uniform mass density 
pi and neutral particles (n) of uniform mass density p n . The 
equation of motion of the electrons can be written as: 



f^(Ve.V)Ve 



E + 



V e x B 



= -Vpc - 



Vn) 



(1) 



where the electron-ion collisions have been neglected since 
the ionized component is of low density. On neglecting the 
electron inertial force, the electric field E is found to be: 



E 



V x B Vp c m c 

V<m(V e - Vn) 

c erie e 



(2) 



This gives us Ohm's law. For 5 = (pi/p n ) <S 1 the ion dy- 
namics can be ignored. The ion force balance then becomes 

Vi x B 



= — Vpi + eni 



E + 



- VinPiiVi - V n ) 



(3) 



where v ln is the ion-neutral collision frequency, and the ion- 
electron collisions have been neglected for the low density 
ionized component. Substituting for E from Eq. (2) we find 
the relative velocity between the ions and the neutrals: 



Vn ~ Vi = 

where 



V(pi+p c ) JxB 



Vinpi 



CVinPi 



J = en c (Vi - V e ) . 

The equation of motion of the neutral fluid is 
-dVn 



(4) 



(5) 



Pn 



dt 



+ (Vn ■ V)V n 



= -Vp n - ^niPn(V n ~ Vi) — 

Wn(K - V e ) , (6) 



where the viscosity of the neutral fluid has been neglected. 
Substituting for V n — V from Eq. (4), and using v m pi = VniPn 
we find: 



Pn 



dVn 

dt 



+ (Vn ■ V)V n 



-Vp + 



JxB 



(7) 



where p — p n + p\ + p c - Observe that the neutral fluid is 
subjected to the Lorentz force as a result of the strong ion- 
neutral coupling due to their collisions. 
Consider Faraday's law of induction: 

OB 

dt 

By substituting for the electric field from Eq. (2) , we get 



= -cV x E 



(8) 



— = V x (V e x B) + riV 2 B , 



(9) 



where the pressure gradient terms have been dropped for the 
incompressible case with constant temperature. Here r\ — 
m e f on c 2 /(47re 2 n e ) is the electrical resistivity predominantly 
due to electron-neutral collisions. Using the construction 



Vc X B = [Vn - (Vn - Vi) - (V - Ve)] X B 



(10) 



and substituting for the relative velocity of the ion and the 
neutral fluid from Eq. (4), Eq. (9) becomes: 



dB 

dt 



= V x 



V n 



.7 JxB 

en c cVinPi 



x B 



+ V X7 2 B (11) 



One can easily identify the Hall term (J/en e ), and the am- 
bipolar diffusion term (J x B) (Chitre & Krishan 2001). The 
Hall term is much larger than the ambipolar term for large 
neutral particle densities or for v\n ^> cJci where io c i is the 
ion cyclotron frequency. In this system the magnetic field is 
not frozen to any of the fluids. Equations (7) and (11) along 
with the mass conservation 

V ■ Vn = (12) 
form the basis of our investigation. 



3 THE ALPHA EFFECT IN 

THREE-COMPONENT MAGNETOFLUID 

The a effect along with its several variants is the key concept 
in the generation of large scale magnetic fields from small 
scale velocity and magnetic fields in the kinematic dynamo 
process (Krause & Radler 1980). The magnetic induction 
equation (11) is written as: 

^ = V x [V E x B] + V X7 2 B , (13) 
where 



V E = K + V H + V Am 
with 

Vk = -± 

en c 

as the Hall velocity and 
JxB 



V A n 



C^inPi 



(14) 



(15) 



(16) 



could be called the ambipolar velocity. Following the stan- 
dard procedure (Krause & Radler 1980) the velocity Ve and 
the magnetic field B are split into their average large scale 
parts and the fluctuating small scale parts as: 



V E = Ve + Ve, 
B = B + B' 



such that 



0, 



0. 



(17) 
(18) 

(19) 



In the kinematic dynamo the magnetic induction equation is 
solved for large and small scale fields. Substituting Eqs. (17) 
and (18) into the induction equation (11), we find, in the first 
order smoothing approximation, 



v , = v ,_£_ + ^ + M 

en c CVinPi CVinPi 

and the mean flow is found to be: 



V E = K - — 

en e 



JxB 



(20) 



(21) 



CVinPi 

The turbulent electromotive force £ is defined as 



£ = Vg X B' and is a function of the mean magnetic induc- 
tion B and mean quantities formed from the fluctuations. 
The fluctuations in turbulence have, generally, a correlation 
in spatial scale L cor and time scale r cor . In a two-scale tur- 
bulence, L CO r <S L and r CO r <§C t, where L and f represent 
the scales of the large scale quantities. Thus the fluctuations 
need to be determined in the immediate vicinity of the point 



Mean-field dynamo in partially ionized plasmas-I 3 



at which the large scale quantity is to be found. This enables 
us to employ Taylor's expansion for B and express the tur- 
bulent electromotive force as (Eq. 5.4 of Krause & Radler 
1980), retaining only the first order spatial derivatives and 
omitting all time derivatives, 



£i= (V^xB') i = a ij .B j +& iik ^ 



(22) 



For a zero mean flow (Ve = 0) , homogeneous, isotropic, 
steady and non-mirror symmetric turbulent velocity field 
Ve> the coefficients ay s and & ; j k become isotropic, and can be 
expressed as: 



dij = aSij ; 



&ijk = /3eijk 

and the electromotive force can be expressed as: 



£ = aB 



IVxB 



(23) 
(24) 

(25) 



The coefficients dij or a and bijk or /3 become functions of 
B for large B as we find here from the contribution of the 
ambipolar diffusion. The quantity a, a pseudoscalar, turns 
out to be the kinetic helicity of the turbulence and is defined 
as: 



« = -j^.(Vx^) 

= a v + an + a Am ■ 
Here, retaining terms quadratic in fluctuations, 



(26) 



(27) 



is the measure of the average kinetic helicity of the neutral 
fluid in the turbulence possessing correlations over time T cor 
and 

2T CO r " 



QH 



3en ( 



j' ■ n' n 



(28) 



represents the contribution of the Hall effect. The coupling 
of the charged components with the neutral fluid is clearly 
manifest through the possible correlation between the cur- 
rent density fluctuations and the vorticity fluctuations of the 
neutral fluid fl' n = V x V n '. The ambipolar term gives rise 
to 



with 



OLA 



OLA B 



2r c< 



-J' x Sl' n 



(29) 



(30) 



as the contribution from the ambipolar diffusion with its es- 
sential nonlinear character manifest through its dependence 
on the average magnetic induction. One also observes that 
the Hall alpha (Eq. 28) requires a component of the fluctu- 
ating current density along the fluctuating vorticity of the 
neutral fluid whereas the ambipolar effect (Eq. 29) thrives on 
the component of the fluctuating current density perpendic- 
ular to the fluctuating vorticity. The turbulent dissipation 
is given by 



/3=^^ 2 =/3 v +/3 H +/3A n 
with 



(31) 
(32) 



as the measure of the average turbulent kinetic energy of the 
neutral fluid in the turbulence possessing correlations over 
time Tcor and 



2r c , 



3en c 



J' ■ Vn 



(33) 



represents the contribution of the Hall effect. The coupling 
of the charged components with the neutral fluid is clearly 
manifest through the possible correlation between the cur- 
rent density fluctuations and the velocity fluctuations of the 
neutral fluid. The ambipolar term furnishes 



/3ah 
/3a 



Pa B 

2r< 



3cp x v in 



J' X Vn 



(34) 
(35) 



with its essential nonlinear character manifest through its 
dependence on the average magnetic induction. One also 
observes that the Hall /?h requires a component of the cur- 
rent density fluctuations along the velocity fluctuations of 
the neutral fluid whereas the ambipolar effect thrives on the 
component of the current density fluctuations perpendicular 
to the velocity fluctuations. We have used rigid or perfectly 
conducting boundary conditions (all surface contributions 
vanish) while determining the averages. Here we consider 
what is known as the a 2 dynamo and take the mean flow 
Ve = 0. This actually determines the relative mean flow 
amongst the three fluids. The dynamo equation reduces to 

dB 



dt 



= V x [aB - /3V x B]+ r]V B 



(36) 



From now onwards we omit the bar on the large scale mag- 
netic field. Assuming one dimensional space dependence, 
we assume the magnetic induction B = (0, B, dA/dx) in 
Cartesian coordinates (a;, y, z). In the corresponding spher- 
ical configuration, one identifies the Cartesian coordinates 
(a;, y, z) with the polar coordinates (8, <f>, r). Thus B and 
dA/dx represent the toroidal and the poloidal parts of the 
field, respectively (Stix 1972). The boundary conditions then 
turn out to be the vanishing of B and A at the endpoints 
of a finite x-interval, say x = and x = irR corresponding 
to the poles of the sphere. It is convenient to put the induc- 
tion equation in a dimensionless form using a normalizing 
magnetic field Bo, a spatial scale R, a time scale R 2 jr\\ and 
writing A = A'R. It begets: 



dB_ 

dt 



and 
OA' 



d 2 A, d 2 B d 
ox ox dx 



B + a 



dA' 
dx 



dA' 

dx 



R/3A 



d_ 

dx 



B + b 



dA' 
dx 



dB 

dx 



(37) 



d 2 A' ( dA n 

= R a B + — — + R aA B + a—— | B + 

dx 2 \ dx 



RfjA [B + b- 



dA'\ d 2 A' 



dx I dx 2 



where 

Ra = 



R, 



olA 



aiR 

a Ay RBo 



(38) 

(39) 
(40) 
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R/3A 



b = 



QAz 
"Ay ' 

/3 Az 



Vi = V + + /3h, 
ai = a v + Qh • 



(41) 
(42) 

(43) 

(44) 
(45) 



Here one observes that since the Hall effect contributes lin- 
early, it can be combined with the standard a v effect. The 
ambipolar effect is nonlinear and appears separately in the 
induction equation. The Hall effect can completely quench 
or enhance the standard a v contribution to the dynamo for 
V n ' = ±J'/(en e ). In the absence of the ambipolar effect one 
recovers the standard a 2 effect with an exponential growth 
rate of the magnetic induction. 

It is instructive to examine the new correlations for the 
case of say Alfvenic turbulence. Now in the weakly ionized 
case the relation between the velocity and the magnetic 
field fluctuations for the Alfven waves is given by = 
±5B'/^/Wpi (Krishan & Varghese 2007) with 8 = pi/p n . 
Substituting these results in the expression for qh we find: 



qh 



^2Ah 



(46) 



where Ah = c/ojpi is the ion inertial scale and w p i = 
\J Annie 2 /m,i is the ion plasma frequency and 



(47) 



is the ratio of the average kinetic helicity and the average 
enstrophy of the neutral fluid turbulence. It is interesting to 
note that the ambipolar a effect vanishes for the Alfvenic 
turbulence. In the next section we discuss the results of the 
numerical solutions of Eqs. (37) and (38) . 



4 DISCUSSION AND CONCLUSION 

We solve the field Eqs. (37) and (38) demonstrating the 
linear and the nonlinear a effect for the initial conditions 
given by Stix (1972): A'(x, 0) = 0, B(x, 0) = sin a;. In the 
absence of the ambipolar contribution (R&A = 0, RpA = 0). 
The Eqs. (37) and (38) become linear with a solution of the 
form exp[i(kx + yt)], where 7 = — k 2 ± kR a and k is the 
dimensionless wavenumber. Then the field grows for k < R a 
i.e., for large values of the effective a at large spatial scale. 

With the inclusion of the ambipolar contributions equa- 
tions become nonlinear. In Fig. 1, we present a case with 
the Hall and the ambipolar contributions. The comparable 
values of the coefficients of the linear and the nonlinear a 
terms with R a = 1.6, -R^a = 1.7, RpA = 0.1, a = 1 and 
b — 0.7 lead to a near constant toroidal field near x = 2, 
fast decaying solution at x = 0.4 and growing solution at 
x — 0.75 beyond t ~ 0.5. The poloidal field, however, grows 
at all values of x. The panels c and d demonstrate the ex- 
pected formation of spatially sharp magnetic structures due 
to the nonlinearity of the ambipolar diffusion (Brandenburg 
& Zweibel 1994). The toroidal field in addition undergoes a 
reversal at x ~ 0.5. Figure 2 shows the dominant effect of the 



ambipolar term with R a = 0.2, -R^a = 3.5, RpA = 1.5, a = 1 
and b = 1. Both the components of the magnetic field, after 
an initial near steady state, grow rather fast and again the 
formation of small spatial scale structures is evident. Thus 
the inclusion of the Hall and the ambipolar effects opens up 
a range of possible profiles of the magnetic field. 

In this first attempt, a framework and some instructive 
examples of the dynamo solutions in a 3-component mag- 
netofiuid have been given. In subsequent work we plan to 
investigate the role of the Hall and the ambipolar terms in 
some realistic situations such as the Solar surface, molecular 
clouds and the accretion disks. In order to deal with these 
systems, the differential rotation in the objects must be in- 
cluded. In the three component system, one would have to 
specify the rotation profile of all the components since the 
system can afford to carry a net current density. The in- 
clusion of the ion-neutral collisions introduces an additional 
time scale with which the turbulence correlation time needs 
to be contrasted. The inclusion of the Hall effect brings in 
the physics at the ion-inertial spatial scale and ion gyration 
time scale. The possibilities are many and varied and should 
be explored in a system specific manner. 
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